Non-Fermi Liquid due to Orbital Fluctuations in Iron Pnictide Superconductors 
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We study the influence of quantum fluctuations on the electron self energy in the normal state 
of iron-pnictide superconductors using a five orbital tight binding model with generalized Hubbard 
on-site interactions. Within a one- loop treatment, we find that an overdamped collective mode 
develops at low frequency in channels associated with quasi-lD dxz and dy^ bands. When the 
critical point for the Ca symmetry broken phase (structural phase transition) is approached, the 
overdamped collective modes soften, and acquire increased spectral weight, resulting in non- Fermi 
liquid behavior at the Fermi surface characterized by the frequency dependence of the imaginary 
part of electron self energy of the form c^;^, < A < 1. We argue that this non- Fermi liquid 
behavior is responsible for the recently observed zero-bias enhancement in the tunneling signal in 
point contact spectroscopy. A key experimental test of this proposal is the absence of the non- Fermi 
liquid behaviour in the hole-doped materials. Our result suggests that quantum criticality plays an 
important role in understanding the normal state properties of iron-pnictide superconductors. 
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Introduction - Whether or not the iron pnictide su- 
perconductors are strongly correlated materials is hotly 
debated. Certainly a clean association of non-Fermi 
liquid behaviour, either experimentally or theoretically, 
with any part of the phase diagram would suffice to set- 
tle this debate. While parallels with the cuprates are 
suggestive [mHI, they have not resulted in a decisive an- 
swer to this problem. In fact, to our knowledge, the pos- 
sibility of non-Fermi liquid behaviour other than Mott 
physics [1] has not been discussed to date. 

A feature common to the parent and undcrdoped 
compounds of the iron pnictide superconductors is the 
structural phase transition (SPT) from tetragonal to or- 
thorhombic symmetry occuring around 150K 7]. For 
most members of the 1111 and 122 families, in-plane 
anisotropy in the resistivity commences near the struc- 
tural transition, and stripe-like antiferromagnetism de- 
velops if the temperature is further lowered. The quasi- 
particle interference in STMj^ also showed the an- 
iostropic electronic states at very low temperature. De- 
spite the controversy of whether the SPT is induced 
by magnetic fluctuations as a result of the onset of 
stripe-like antiferromagnetism l^- 14 1 or if orbital order- 
ing in quasi-lD dxz and dyz bands isl- 18 1 is the effi- 
cient cause, the phase below the SPT breaks C4 sym- 
metry, and quantum fluctuations associated with this 
phase are nematic in character. A recent measurement 
of photoexcited quasiparticle relaxation[l^ reveals the 
existence of strong nematic fluctuations up to 200K, 
well above the SPT temperature. Moreover, in electron- 
doped Ba(Fei_a;Co2;)2As2 (Bal22) an unexpected en- 
hancement of the zero-bias signal [2o[[2l| in the tunneling 
signal measured in point contact spectroscopy has been 
observed at an onset temperature higher than the SPT 
temperature. The excess conductance appears at tem- 
perature around 175K, increasing in magnitude through 
the structural, antiferromagnetic, and, in materials ex- 



hibiting superconductivity, through the superconducting 
transitions. It is not seen in overdoped Bal22. It is im- 
portant to study how these strong orbital (nematic) fluc- 
tuations affect the physical properties in both the normal 
and orthorhombic states of the iron pnictide supercon- 
ductors. 

In this Letter, we develop a microscopic theory for the 
orbital fluctuations and show that they give rise to non- 
Fermi liquid behaviour. In particular, we flnd a branch of 
overdamped collective modes in the scattering channels 
associated with quasi-lD dxz and dyz bands in the nor- 
mal state at a temperature higher than the SPT. In the 
vicinity of the SPT critical point, these overdamped col- 
lective modes dominate the low energy physics, resulting 
in a strong modiflcation of the electron self energy and a 
breakdown of Fermi-liquid theory even in the symmetric 
normal state. 

RPA theory for electron self energy - While there is no 
guarantee that non-Fermi liquid behaviour can be cap- 
tured by a diagramatic approach of the type employed 
here, a recent non-perturbvative calculation of the elec- 
tron self-energy in the presence of nematic fluctuations 
was found to be in full agreement with a perturbative 
treatment [i^l . As the same physics is relevant here, we 
contend that a perturbative treatment is sufficient to es- 
tablish our key claim of non- Fermi liquid behaviour from 
orbital fluctuations. Our starting point is the Hamilto- 
nian H = Ht + Hj, where Ht — is the five-band 

tight binding model proposed in Ref. 23] which can re- 
produce correctly the Fermi surfaces of hole 01,0:2 and 
electron pockets /3i , (52 in the unfolded Brillouin zone. 
The interaction terms are given by Hj 
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FIG. 1: Spectral functions for effective interactions in intra- 
orbital channels {raa;aa{q,ijj)) at qi = (0.047r, 0.047r) for (a) 
U=1.3 and (b) U=2.0. rii;ii(g, w) = T22;22{q,oj) for q along 
the diagonal direction as expected, and an overdamped collec- 
tive mode appearing at low energy in rii;ii(if, aj)(r22;22(g, t^)) 
can be observed. 



where U' = U — 2J and J is the Hund coupling, and 
a (b) refers to the orbital index, 1 — xz, 2 = yz, 



xy, 
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We introduce 



a unitary transformation JJj^ such that (Uj^ k ~ 
diag. \_Ej: j^, • • • , Ej: , and it is straightforward to obtain 
the non-interacting response functions. 
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in the symmetric normal phase, where xft-cdi'!^'''^"-) is a 
25 X 25 matrix. The convention on the indices is that in 
Ref. [i^l. Adopting the interaction kernels for the spin- 
spin (V_f) and density-density {V'^) fluctuations derived 
in Ref . |24| . we obtain the electron self energy at the 



FIG. 2: Spectral functions of effective interactions 
rii;ii((fi,tj) at (fi = (0.047r, 0.047r) for different U. The over- 
damped collective mode at low energy acquires increased spec- 
tral weight as the critical point {Uc ~ 2.feV at the tempera- 
ture we are considering fcaT — 0.02eV) is approached. 



random-phase approximation level, 

Kl'^'^'ik, iPn) = ^ E E iPn-lkrr.)Gl{k-q, tk„ 

q ikm 

(3) 

where ikm,,ipn are Matsubara frequencies for fermions, 
&{k, ipn) = \}Pn — Ht + lj\ is the bare Green function, 
and T ab;cd{<l, iPn — ikm) is the effective interaction, 

t{q,ipn - ikm) = i {3 [l - 1^^X0(9, «Pn - «fcm)] ^V" 
- [l + V''X0{q,Vn'~lkm)r'V'} (4) 

within one-loop. Throughout this paper, the tempera- 
ture is set to be ksT — 0.02eV at which the system is in 
a normal state without any symmetry breaking. 

Consider first the spectral functions for the effective 
interactions —lmT{q,uj + irf) displayed in Fig. [TJ As is 
evident, the spectral functions for the intra-orbital ef- 
fective interactions —ImT aa\aa(<li,'-^) dominate the elec- 
tron self energy. Also of interest are the spectral func- 
tions for momenta q along the diagonal direction so that 
riiiiiCif, w) = r22;22(9, w). It can be seen in Fig. [T]that 
the spectral functions at low frequency are dominated by 
an overdamped collective mode in Fn;!! (F22;22). When 
U is tuned to approach the critical point [Uc ~ 2.1 at 
this temperature), this mode gains more spectral weight 
and moves to even lower energy as shown in Fig. [51 

These overdamped modes, emergent at low frequency 
and small q, resemble the collective modes observed in 
the quadrupole density spectral function 22 , H, [2^ , that 
is, the spectral function related to the interactions in the 
d-wave channel in a quantum nematic Fermi fluid. As 
shown in Ref. f23l, in a system containing quasi-lD dxz 
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and dyz bands, band hybridizations enhance significantly 
the strength of the interaction in the cZ-wave channel As 
a result, the nematic order in such multiorbital systems 
is completely equivalent to orbital ordering in quasi-lD 
bands, and the spectral functions due to quantum fluctu- 
ations associated with the quasi-lD bands naturally ac- 
quire the same properties of quadrupole density spectral 
function discussed in the context of the quantum nematic 
fluid mentioned above. It has been shown 22, 25, 2^ that 
these overdamped collective mode could lead to a non- 
Fermi liquid near the critical region and also in the ne- 
matic phase. The reason is that in the vicinity of the ne- 
matic critical point, these overdamped collective modes 
become soft. Electrons scatter strongly with these soft 
overdamped collective modes, which modifies the elec- 
tron self energy away from the Fermi liquid behavior 
in the vicinity of the nematic critical point. In the ne- 
matic phase, these overdamped collective modes evolve 
into Goldstonc modes [2^ but remain overdamped and 
dominate the low energy physics. Consequently, the non- 
Fermi liquid persists in the nematic phase. 

It is intriguing to check whether the same physics dis- 
cussed above occurs in the iron pnictide superconductors 
since the SPT signals a transition from the symmetric 
normal phase to a phase breaking C4 symmetry. We per- 
formed a numerical evaluation of Eq. [3l To compute the 
self energy of the retarded Green function of a single- 
particle state on the Fermi surface, we need to do one 
more transformation and also an analytical continuation 
on Eq. |3]to obtain 
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which is the self energy of the electron with momen- 
tum kp on the Fermi surface sheet a. Because we study 
the normal state at finite temperature, T,^^'^{kp, uj + iij) 
contains contributions from both thermal and quantum 
fluctuations. At low frequency, we can generally write 
-ImS^"„"'^(fcF, w + iri) « a(T) + b{T)oj^ where a{T), b{T) 
are temperature-dependent constants. A Fermi liquid 
yields A = 2, and A < 1 signals a non-Fermi liquid. 

To see the frequency dependence more clearly, in Fig. 
[3]wc plot the normalized imaginary part of the self energy 
defined as 



ImSS°'-(fcF,w) 



(lmE^'^„"'^(fcf,c^) 
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(6) 

A crossover from Fermi liquid (A = 2) to non- Fermi liq- 
uid (A < 1) behaviour can be seen as U is increased to 
approach the critical point. In the critical region, non- 
Fermi liquid behavior exists in a large part of the Fermi 
surface with strong angular dependence of ImS as ex- 
pected due to the critical fluctuations near orbital order- 
ing (now termed nematicitv) psf . As the temperature is 
lowered, the critical point Uc shifts to a lower value but 
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FIG. 3: (a)Normalized self-energy of electron with momen- 
tum kp = (0.127r, 0.127r) on the hole Fermi pocket qi for 
different U. A crossover from Fermi liquid (A = 2) to 
non-Fermi liquid (A < 1) can be seen as U increases from 
U = 0.5 to critical point U = 2.1. (b) Normalized self- 
energy of electron in the critical region (U=2.1) with mo- 
menta kpi = (0.127r, 0.127r), kpi = (0.27r,0) on hole Fermi 
pocket Qi, and fc^i ~ (O.SStt, O.IGtt) on electron Fermi pocket 
/3i. 



the non-Fermi liquid behavior remains robust near the 
critical point. This strongly suggest that this non-Fermi 
liquid behavior should be observable in iron-pnictide su- 
perconductors at a temperature above the SPT and will 
persist even below the SPT temperature. A sketch of the 
region in the phase diagram exhibiting the non- Fermi liq- 
uid behavior is plotted in Fig. HI 

Experimental Consequences - A direct cconsequence 
of non-Fermi liquid behavior is the temperature depen- 
dence of the resistivity. It has been pointed out [2^, by 
investigating various iron-pnictide superconductors, that 
a strong deviation from the Fermi liquid behavior of 
resistivity above the SPT temperature would occur if a 
large anisotropy in in-plane resistivity exists below the 
SPT temperature. This provides direct evidence for our 
claim that non-Fermi liquid behaviour is due to orbital 
fluctuations. 

What about the zero-bias anomaly seen in point- 
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FIG. 4: (Color online) Sketch of regions in the phase diagram 
showing non- Fermi liquid behavior (colored area). Solid line 
marks the phase boundary of the transition from symmetric 
to C4 symmetry breaking phases. 



modes emergent at low frequency in channels associated 
with quasi-lD dxz and dyz bands, and we have shown 
that these modes become dominant at low energies near 
the critical point for C4 symmetry broken phase, lead- 
ing to non-Fermi liquid behavior. Our theory indicates 
that quantum criticality could play an important role 
in understanding the normal state properties of the iron 
pnictide superconductors. 
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contact experiments [20|, |2l| on electron-doped 
Ba(Fei_a;Coa;)2As2? Intriguingly, this zero-bias en- 
hancement starts to appear at temperatures roughly 
30K higher than the SPT temperature and remains ro- 
bust well below the SPT temperature. This observation 
is also consistent with our theory. It has been shown 
by Lawler et. al.^^ that the single particle density of 
states has the form of 



N*{lu) ^ N*{0) + Blu"^/^ In a; + 



(7) 



in the nematic critical region and also in the nematic 
phase. In fact, N*{oj) obtains extra contributions due 
to the non-Fermi liquid self-energy, giving rise to a cusp 
at zero frequency and a subsequent decrease as the fre- 
quency increases. This provides a direct explanation 
for the zero-bias enhancement observed in point contact 
spectroscopy since the conductance dl/dV roughly mea- 
sures the the single particle density of states for small 
frequency. Moreover, since the form of the single particle 
density of states is the same up to some mild modifica- 
tions in the vicinity of the critical point and also in the 
C4 symmetry broken phase, the zero-bias enhancement 
should have a smooth crossover as the SPT temperature 
is crossed, which in fact has been noticed in the quantum 
point-contact measurement [iS [U. We predict that for 
hole-underdoped Bai_a;Ka;Fe2As2 which does not have 
an in-plane resistivity anisotropyfl^, the zero-bias en- 
hancement should be either not existant or much weaker 
than that in electron-doped Ba(Fei_2:Co2:)2As2. 

Conclusion - In this Letter we have presented a the- 
ory of non-Fermi liquid in a five-band model with gen- 
eralized Hubbard on-site interactions for iron pnictide 
superconductors. At the level of the random-phase ap- 
proximation, we found a branch of overdamped collective 
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